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ABSTRACT 

During its 2005 outburst, GRO J1655-40 was observed twice with the Chandra High Energy Trans- 
mission Grating Spectrometer; the second o bservation revealed a spectrum rich with ionized absorption 
lines from elements ranging from O to Ni (IMiller et al.M2006al [200l IKallman et all [2001 . indicative 
of an outflow too dense and too ionized to be driven by radiation or thermal pressure. To date, this 
spectrum is the only definitive evidence of an ionized wind driven off the accretion disk by magnetic 
processes in a black hole X-ray binary. Here we present our detailed spectral analysis of the first Chan- 
dra observation, nearly three weeks earlier, in which the only signature of the wind is the FeXXVi 
absorption line. Comparing the broadband X-ray spectra via photoionization models, we argue that 
the differences in the Chandra spectra cannot possibly be explained by the changes in the ionizing 
spectrum, which implies that the properties of the wind cannot be constant throughout the outburst. 
We explore physical scenarios for the changes in the wind, which we suggest may begin as a hybrid 
MHD/thermal wind, but evolves over the course of weeks into two distinct outflows with different 
properties. We discuss the implications of our results for the links between the state of the accretion 
flow and the presence of transient disk winds. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: individual (GRO J1655- 
40) — X-rays: binaries — stars: winds, outflows 



1. INTRODUCTION 

The last two decades have seen the discovery of 
a multitude of highly-ionized absorbers and winds in 
the X-ray spectra of accr eting black hol e and neutron 
star X-ray binaries (e.g . lEbisawal Il997t IKotani et al.l 
19971 iBrandt fc Schulzl 120001: IKotani et al J l2000allb : 
Lee et al.ll2002t ISchulz fc Brandtlte002l: lUeda et al 1 1 2004 : 
20041 l2006aUbl. 120081: INeilsen fc Led I2OO9I: 
2009fc iSTim et all 120101: IRevnolds fc Miller! 
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Kin g et alJl20T2ir In many of these cases, although the 



kinetic power in these outflows is orders of magnitude be- 
low the accretion power, the mass loss rate in the wind 
can be much greater than the accretion r ate at the inner 
edge of the disk (e.g. INeilsen et alJl2011f ). The ubiquity 
of disk winds in X-ray binaries and their extreme mass 
loss rates suggests that these outflows play a crucial role 
in the physics of accretion. 

Presently, accretion disk winds are understood to be 
driven by some combination of radiation, thermal, and 
magnetic pressure, which operate in different regimes 
of density and ionization (iBegelman et al.l H983: Proga 
[20001: IProga fc Kallmanl I2002t IMiller et al.ll2006afl . Yet 
determining which of these processes dominates the 
physics of an individual outflow typically presents a sig- 
nificant challenge. Most often, indirect density estimates 
must come from geometrical arguments that constrain 
the distance R of the wind from the black hole (e.g. 
iLee et al.|[2002l: INeilsen et al.ll201lL [20121 ) . In these cases, 
the density n can be found usin g the observed ion ization 
parameter £ and luminosity L ([Tarter et al.1 11969): 



L 
~nR?' 



(1) 



In rare cases it may be possible to measure the den- 
sity directly from density-sensitive atomic lines, a tech- 
nique that was most notably applied in the d iscovery of a 
magne tically-dri ven wi nd in GRO J1655-40 ([Miller et al.l 
l2006al hereafter [M06). When Chandra observed this 
microquasar with the High Energy Transmission Grat- 
ing Spectrometer (HETGS; iCanizares et al.ll2005T) in the 
middle of its 2005 outburst, the high-resolution X-ray 
spectrum (ObsID 5461) revealed a rich series of absorp- 
tion lines from a dense (n > 10 14 cm" 3 ), highly- ionized 
(£ ~ 10 4 ) wind. Detaile d follow-up spect ral analy- 
sis and theoretical studies ([Miller et al.1 120081 hereafter 
MOH IKallman et all [20091 . hereafter |K0l iLuketic et al.l 
20101 ) have repeatedly confirmed the original conclusion 
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of lMOft this outflow was driven primarily via magnetic 
processes, rather than radiation or thermal pressure. 

However, a HETGS observation (ObsID 5460) of the 
same outburst made 20 days earlier reveals an outflow 
that is empirically quite different (see Figure [1]). Instead 
of the host of lines, this observation merel y cont ains an 
FeXXVi absorption line noted in passing by IM08I . In this 
paper, we ask: where did the other lines go? During the 
first observation, are we simply observing a much more 
ionize d (and thus tra nspar ent) version of the wind stud- 
ied bv lM06llM08l . and lK09l ? Or do these two observations 
actually probe a transient wind that evolves physically 
and geometrically throughout the black hole outburst? 

We describe the observations and data analysis in Sec- 
tion [2J In Section [31 we present the Chandra spectra and 
use broadband X-ray spectra from the Rossi X-ray Tim- 
ing Explorer's (jJahoda et al.|[l996l) Proportional Counter 
Array (RXTE PCA) to study the ionizing continuum. 
We also build photoionization models to demonstrate 
how the wind responds to changes in the incident X-ray 
spectrum. In Section [4[ we discuss the implications for 
the evolution of wind physics during the 2005 outburst. 
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Fig. 1— Chandra HETGS spectra of GRO J1655-40 during its 
2005 outburst. Shown in black is the spectrum of the first obser- 
vation, when the broadband X-ray continuum was hard and very 
few absorption lines were present; shown in blue is the spectrum 
of the second observation, which was executed 20 days later and 
recorded a rich series of lines from the accretion disk wind. In this 
paper, we address the physics driving these significant changes in 
the disk wind. 

We conclude in Section [5j 

2. OBSERVATIONS AND DATA REDUCTION 

GRO J 1655-40 was observed with the Chandra 
HETGS on 2005 March 12 (20:42:53 UT; ObsID 5460) 
and 2005 April 1 (12:40:41 UT; ObsID 5461) for 34.3 and 
62.15 ks, respectively. The data were taken in Graded 
Continuous Clocking mode, which has a time resolution 
of 2.85 ms and is designed to limit both photon pileup 
and telemetry saturation. 

We reduce and barycenter-correct the Chandra data 
with standard tools from the CIAO analysis suit, version 
4.3. We use the order-sorting routine to remove the read- 
out streak on the chip S4, since the dedicated tool (de- 
streak) can also remove source counts when the count 
rate is high. After reprocessing, we extract High-Energy 
Grating (HEG) spectra and make grating response files. 

During the Chandra observations, RXTE made 
pointed observations of GRO J1655-40 with good ex- 
posure times of 1.82 and 12.08 ks, respectively. In this 
paper, we analyze data from the PCA, which covers 
the 2-60 keV band, in order to explore the role of the 
photoionizing continuum in driving the evolution of the 
wind. For this spectral analysis, we extract Standard-2 
129-channel spectra, and restrict our attention to the 3- 
45 keV spectrum from PCU2. We assume 0.6% system- 
atic errors in ea ch energy bin. All our spectral fittin g 
is done in ISIS (IHouck fc Denicolal [20001 iHouckl l200l . 
Following iMOl we assume a neutral hydrogen column 
density A H = 7.4 x 10 21 cm -2 and a distance of 3.2 kpc 
(|Orosz fc Bailvnlfl99l . 

3. X-RAY SPECTROSCOPY 

The Chandra HEG spectra are shown in Figure [TJ 
where it is abundantly clear that there are very signif- 
icant differences between the absorber during the first 
(black) and second (blue) observations. In this section, 
we characterize the high-resolution and broadband spec- 
tra in an effort to determine whether photoionization 
alone can explain explain the differences in the lines. Our 
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Fig. 2.— RXTE PCA spectra of GRO J1655-40 during its 2005 
outburst. Shown in black is the spectrum of the first observation, 
when the broadband X-ray continuum was hard and very few ab- 
sorption lines were present; shown in blue is the spectrum of the 
second observation, which occurred 20 days later and exhibits a 
significantly softer X-ray continuum. In Section [3] we argue that 
the changes in photoionizing flux are insufficient to explain the 
changes in the lines (Figure [TJ . 

method is as follows: first we model the HEG/PCA spec- 
tra, then use xstar to calculate the expected atomic 
level populations for a gas slab illuminated by our ob- 
served ionizing continuum. Finally, we compare the pre- 
dicted and observed absorption spectra. We focus on the 
first observation, which has not been studied in detail. 

3.1. Preliminary Modeling 

To start, we fit the HEG spectra phenomenologically 
with an absorbed polynomial spectrum and multiplica- 
tive Gaussian absorption lines (multiplicative absorption 
models are accurate even if the lines are moderately opti- 
cally thick; Hdi^irSIlQil). The primary goal of this 
phenomenological modeling is not to catalog the lines 
or investigate their properties in detail (which has been 
done several times for the second observation, e.g. IM081 
IK09T ) . but to incorporate the ionized lines into models of 
the broadband X-ray spectrum as seen by RXTE. Thus 
we devote relatively little space in this section to a de- 
tailed discussion of the lines themselves. 

However, a direct comparison of the FeXXVl Lya lines 
in the two observations is informative. In both spectra, 
we achieve tighter constraints on the parameters of these 
lines by fitting the Lya and Ly/3 lines simultaneously 
(i.e. with the same ion column density, velocity, and line 
width), since the optical depth in a line from any given 
ion is proportional to the line oscillator strength times 
the ion column density divided by the line width. Thus, 
even though Ly/3 is not observed during observation 1, 
its absence places an upper limit on the column density 
of Fexxvi. In the first observation, we find a column 
density Ap CX xvi = 6-1 x 10 17 cm -2 and a blueshift of v — 
— 470±230 km s _1 (the errors are 90% confidence limits). 
The line width is a = 13^ eV, and the equivalent width 
is Wo ~ 24 eV. In the second observation, the FeXXVl 
has a similar column density, Apcxxvi = (7.6 ±0.9) x 10 17 
cm~ 2 , but a higher velocity, width, and equivalent width 
(v = -1350t^ km s" 1 , a = 21 ± 4 eV, and W Q w 31 
eV). In the second observation, at least 33 additional 
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strong features are present (in the 2.5-9 keV range) that 
are not visible in the first dataset. 

In short, we see that the Fexxvi column densities are 
similar, but all the other ion column densities increase 
dramatically from the first observation to the second, 20 
days later. The simplest explanation for this fact is that 
the wind is more ionized during the earlier pointing. To 
test this hypothesis in more detail, we turn our atten- 
tion to the broadband PCA spectra, which are shown in 
Figure [21 From this figure, the changes in the ionizing 
continua between the two observations are immediately 
apparent: during the first observation, there were many 
more photons capable of ionizing hydrogen- and helium- 
like iron (i.e. with E > 9 keV). On the other hand, other 
elements have very small photoionization cross-sections 
for E > 9 keV. So if the geometry or bulk properties of 
the wind were constant in time, we would expect there 
to be relatively little difference between the two obser- 
vations in terms of absorption lines from ions other than 
iron. Since we observe highly-significant changes in lines 
other than iron, we can already conclude that excess ion- 
izing luminosity alone is definitely not responsible for the 
absence of lines in the first observation. 

Our best fit model for the broadband continuum 
in the first observation (tbabs* (simpl ®ezdiskbb 
+egauss) *nmgauss*edges) consist s of a n absorbed 
(tbabs; Wil ms. Allen, k McCravl 120001) hot disk 
(ezdiskbb; iZimmerman et aLl |2005|) convolve d through 
a scattering kernel (simpl: iSteiner et al.ll2009f ). a Gaus- 
sian emission line at 6.4 keV, along with th e Fexxvi 
absorption line (nmgauss ; iNeilsen et all I2012T ) detected 
in the HEG spectrum, simpl takes a seed spectrum 
and scatters a fraction fsc of the source photons into 
a power law, approximating the high-temperature, 
low optical dep t h reg ime of Comptonization. Follow- 
ing iSaito et"aLl ()2006l ). we allow for possible ionized 
absorption edges at 7.7, 8.83, 9.28, and 10.76 keV, 
but none of these are detected at 90% confidence- 
Assu ming a disk inclination of 67° ()Orosz k Bailvnl 
119971 ) and a color correction factor of 1.7, the disk 
normalization N = 390 ± 20 implies an inner disk 
radius R in = 29.lt°; 7 km, or roughly 2 gravitational 
radii, and the disk temperature is T = 0.973 ± 0.009 
keV. We find that (5.6 ± 0.3)% of the disk photons are 
scattered into the power law, which has photon index 
T = 2.11 ±0.04. The implied unabsorbed luminosity 
(10" 3 -10 3 keV) is L = 6.4 x 10 37 ergs s" 1 , roughly 7% 
of the Eddington limit. For the second observation, we 
include a highecut term to account for curvature in the 
spectrum; we can also replace simpliKiezdiskbb with 
ezdi skbb+nthcomp (Zdzia rski et al.l 119961 : iZycki et al.l 
[19991) . For our purposes here, only the luminosity of the 
second observation is relevant: L pa 5.7 x 10 37 ergs s" 1 . 

3.2. Photoionization Modeling 

The immediate goal of our analysis is to determine 
why there are so few lines in the first observation. Are 
the changes in the luminosity and spectral shape enough 
to wipe out the lines? In Section 13.11 we argued that 
the answer is no, and that changes in the wind geometry 
or density are required. In this section, we make this 
argument quantitative with photoionization models. 

The method itself is straightforward. We use XSTAR 



to calculate atomic data for an illuminated slab of gas. 
We define this slab of gas to match the properties of 
the w ind seen in the second observation, as modeled by 
lK09l and we illuminate it with the broadband X-ray 
spectrum of the first observation. In this way, we can 
test whether the same wind, ionized by a harder X-ray 
spectrum, should have been visible in the first d ataset . 
For reference, we base our models on Model 6 of KOfi 
a density of n — 10 15 cm~ 3 , an equivalent column den- 
sity 2Vh = 10 24 cm~ 2 , a 37% partial covering fraction, a 
blueshift of v — 375 km s , and a turbulent line width 
of wturb = 200 km s . We increase the ionization param- 
eter from £ = 10 4 to ^ = 10 4 05 to account for the higher 
luminosity during the first observation. Finally, we use 
the elemental abundances reported in Figure 21 oflKOl 
We incorporate the resulting atomic data into new 
models of the HETGS observations using the analytic fit 
function warmabs. Using parameters equivalent to those 
measured bv lK09t but with an ionization balance appro- 
priate for the first dataset, we find that there are a num- 
ber of lines that would have been detected at high signif- 
icance if the bulk pro pertie s of the wind were the same 
as those measured by IK09I (see Figure [3]) . Three of the 
strongest features predicted by this test model but not 
observed in our data are Sixiv at 2.0, Sxv at 2.6 keV, 
and Ar xvin at 3.2 keV. To give an idea of how strong 
the features are, we note that without introducing any 
additional curvature into the spectrum (only lines), this 
test model increases \ 2 by over 1800 relative to a model 
with no wind at all. This result confirms our earlier con- 
clusion: changes in the spectral shape and luminosity 
cannot be solely responsible for the absence of spectral 
lines in the first observation. 

4. DISCUSSION 

The results of our photoionization models in Section 
13.21 are unequivocal: if the bulk properties of the wind 
were constant in time during the 2005 outburst of GRO 
J1655-40, a large number of strong lines should have 
been visible in the first dataset. A constant wind, sim- 
ply over-ionized in the hard state, is ruled out. Fur- 
thermore, since the observations occur at orbital phases 
0.5 and 0.0 (where t he disk would be eclipsed at 0.75; 
IQrosz k Bailvn|[l997l) . orbital or companion star effects 
seem unlikely to be important. These results have broad 
implications for the role of winds in outburst (e.g. in 
quenching jets; INeilsen k Led [2009). In this period of 
three weeks, as the source moved out of the hard, jet- 
producing state into a spectrally-soft state, the structure, 
density, or geometry of the accretion disk wind must have 
changed significantly. But what else can be said about 
the details of this evolving wind, given that so few lines 
are visible in the first dataset? 

Most of what we know for certain comes from observa- 
tion 2, which shows not only clear indicators of an MHD 
wind, but a lso compelling eviden ce for a second wind 
component: iKallman et al.l (|2009l ) found that while the 
typical blueshift of absorption lines in the spectrum is 
~ 400 km s -1 , the strong lines of Fexxvi and Nixxvn 
are blueshifted by ~ 1300 km s _1 . See their Figure 15 for 

2 We compare our results to Kallman et al. (2009) rather than 
the original models of Miller et al. ( 2006a, 2008) because the newer 
models (at a higher density) provide a superior fit to the spectrum. 
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TABLE l 

Photoionization Models of GRO J1655-40 



Energy (keV) 

Fig. 3. — Chandra HETGS spectrum of the first observation 
of GRO J16 55— 40 (black) with m odel (red) based on the wind 
described by Kallman ct al. (2009). If the wind were the same 
during both observations, a number of strong lines would have 
been seen during the first observation, despite its harder spectrum. 
That these lines are not observed indicates that the wind evolved 
significantly throughout the outburst. 

a remarkable demonstration of this result, which they in- 
terpreted as indicative of a second, extremely highly ion- 
ized component. In fact, similar (although less extreme) 
differences in line dynamics have also been seen in other 
X-ray binaries. For example, studying an equally rich 
absorp tion line spectrum in GRS 1915+105. lUeda et al.l 
( 2009) found that while most lines were narrow and mod- 
erately blueshifted (er ~ 70 km s^ 1 , v ~ 150 km s -1 ), 
the FeXXVl line was broader and faster (er ~ 200 km 
s _1 , v ~ 500 km s _1 ). Thus we are convinced that there 
is a two-component outflow in the second observation of 
GRO J1655-40. Since the magnetically-driven wind is a 
dense outflow at small radii, we suppose that the second 
component is a relatively low-density (n < 10 12 cm~ 3 ) 
outflow at large radii (R > 10 11 cm), i.e. as produced by 
thermal driving. 

But how does this color our understanding of observa- 
tion 1? Any explanation must hinge on the Fexxvi line, 
which was present in both pointings. We focus on the 
bulk properties of the wind (i.e. density, location, and 
so on) in the hopes of determining the contributions of 
thermal and magnetic processes to the outflows in GRO 
J1655-40 (radiation pressure is ruled out for the same 
reason as in observation 2: it is ineffective at such high 
ionization parameters). So in short, it comes down to 
whether we believe the iron absorption line in the first 
spectrum should be associated with the slow, dense MHD 
component or the faster, more ionized thermal compo- 
nent in the second spectrum. Thus the interpretation of 
this line can be roughly separated into three scenarios, 
in which the wind in observation 1 is a purely thermal 
wind (Section 14. ip . a magnetically-driven wind (already 
discussed in Section. ET2"]l . or a hybrid thermal/MHD wind 
(Section S3). 

While there is a very large body of literature on the 
theoretical aspects o f the origin and properties of accre 
tion disk winds (e.g. iBegelman et al.l 1983; Shield s et al 
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Note. — n is the electron density input for XSTAR in cm -3 ; N H 
and £ are the best-fit column density and ionization parameter 
for observation 1 in cm~ 2 and ergs cm s~ 1 ; jiar and -Rak are 
the density (cm -3 ) and radius (cm) implied by the corresponding 
value of Nn if AR is fixed for the wind (Section I4.2.2B ; n^ R ^ 
and Rar/r arc the density and radius implied by Nn if AR./R. is 
fixed (Section ETOl . 

formation about which wind processes are expected 
to be most important at any given phase of a black 
hole outburst (cf. t he unified model o f disk-jet cou- 
pling in outburst; [Fender et al.l [2004). For exam- 
ple, are there links between the accretion state and 
the relative importance of radiatively-, th ermally-, and 
magnetically-driven win ds? Observations (jNeilsen et al.l 
120121 iPonti et al.l 12012( 1 suggest that varying illumina- 
tion and shadowing of the outer disk may be important in 
producing transient winds, but it is unclear how the den- 
sity, launch radius, and geometry of these winds (MHD 
or thermal) should evolve over time. In the following 
subsections, we explore the implications of a variety of 
simple assumptions about the behavior of the wind in 
GRO J1655-40. By no means do these constitute an ex- 
haustive list, but they should be representative of the 
possibilities. 

4.1. Case 1: A Thermal Wind in Observation 1 

Here, in light of the fact that FeXXVl is present in 
both observations of GRO J1655-40, and may arise in 
a thermal wind in observation 2, we consider the pos- 
sibility that this line also arises in a thermal wind dur- 
ing observation 1. In other words, we suppose that the 
thermal wind is persistent, while the MHD outflow is 
transient. This scenario is illustrated in Figure 0] With- 
out decent estimates of the density or ionization pa- 
rameter of the thermally-driven wind, it is very diffi- 
cult to estimate its location. Thermally-driven winds 
can be launched anywhere outsi de 0.1 Compton radii 

Begelman etHI 119831: 



(Rc = 9.8 x IO^ MbhTtc 1 cm 



iWoods et all 119961: iRahoui et all 



20lb|), whereM B H is 



the black hole mass in solar masses and the Compton 
temperature in Kelvins is given by: 



Tic 



1 / °° hvL v dv 
4fc B / °° ' L v dv 



(2) 



1986 ) [Woods et al.llT996t |Progall2000t iProga fc Kallman 
2002t iLuketic et al.l 120101 ). there is relatively little in- 



where h is Planck's constant, v is the frequency, and L v is 
the monochromatic luminosity. For our PCA spectra of 
observations 1 and 2, we find Tic ~ 8.6 x 10 7 K and 8.4 x 
10 6 K, respectively, which imply Compton radii Rq ~ 8x 
10 10 cm and 8 x 10 11 cm. The Compton temperature of 
the first observation is much higher because the spectrum 
is much harder; accordingly, Compton heating is more 
efficient and a thermal wind can arise at smaller radii. 
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Fig. 4. — Cartoon of the scenario presented in Section 14,11 in 
which the wind in observation 1 is thermally-driven (top), while 
observation 2 (bottom) features a two-component wind, i.e. MHD 
and thermal. For the winds, darker shading implies higher density, 
and the different orientations of the two components reflects our 
uncertainty about the vertical component of the wind velocity. See 
text for details. 

Importantly, in both cases, 0.1i?c is inside the disk, so 
that a Compton-heated wind is plausible at this radius. 

The wind speed provides a different perspective. If 
we suppose that the wind speed is comparable to the 
local escape speed in the disk, then the thermal wind 
originates at 10 12 cm and 10 11 cm in the first and sec- 
ond observations, respectively. It is puzzling that the 
implied distance between the wind and the black hole 
in observation 1 is comparable to the binary separation. 
However, the wind velocity likely has a significant com- 
ponent perpendicular to the line of sight, so that the 
Doppler shift represents only a fraction of the true wind 
speed. But even if this is the case, the acceleration of the 
wind still requires some explanation. Why does the wind 
accelerate as the spectrum softens and the luminosity de- 
creases? One possibility tha t has gained some ground 
recently (jNeilsen et al.l (2012t iPonti et all l2012t) is that 
durin g harder st ates, the geometrically-thick inner disk 
(e.g. [Meier 2001) shadows some of the outer disk. This 
shadowing interpretation is potentially testable via opti- 
cal/infrared measures of the irradiation of the companion 
star or outer disk (a paper on multiwavelength observa- 
tions of this outburst is forthcoming). In this framework, 
as the scale height of the inner disk decreases during the 
state transition in GRO J1655-40, progressively smaller 
radii in the disk are irradiated, leading to a faster wind. 
A different irradiation geometry could also influence the 
direction of the wind velocity, e.g. driving it more into 
the line of sight in observation 2. 

In summary, taken at face value, the estimates of R 
from the Compton radius and the escape velocity are not 
fully consistent. However, it may be possible to reconcile 
them by considering the vertical component of the wind 
velocity (which may be significant). Without a clear and 
unambiguous way to estimate the location or the density 
of this thermal wind, there is little more we can say about 
it. Nevertheless, we feel our estimates here are still mer- 
ited, as there is com pelling evidence for a two-component 
wind 20 days later (jKallman et al.ll2009f ). 

4.2. Case 2: An Evolving Hybrid Wind 

On the other hand, given that observation 2 provides 
the only direct evidence for a magnetically-driven wind 
in an X-ray binary, it is also worth considering the pos- 
sibility that both wind components (MHD and thermal) 
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Fig. 5. — Cartoon of the scenario presented in Section 14.21 in 
which observation 1 (top) probes a hybrid wind driven by a combi- 
nation of MHD and Compton heating. In this interpretation, the 
hybrid wind evolves over the course of weeks into two distinct out- 
flows: a thermal wind at large radii and a MHD wind at small radii 
(bottom). See text for details, and see Figure [4] for notes about 
shading and orientation in these cartoons. 

are present in both observations of GRO J1655-40. In 
this scenario, illustrated in Figure [31 a precursor to the 
MHD wind makes a weak contribution to the FeXXVl 
line in observation 1 (which is dominated by the thermal 
wind), but in the subs equen t weeks evolves into the dense 
outflow discovered bv lM06l There is some circumstantial 
evidence to suggest that this scenario is plausible. 

In particular, observations with XMM-Newton on 2005 
March 18 and March 27 (i.e. six days after Chandra ob- 
servation 1 and five days before Chandra observation 
2, respectively) reveal the presence of an absorber with 
an ionization parameter that decreases rapidly over time 
(Diaz Trig o et all l2007f ) . Possible absorption lin es were 
also detected with Swift (jBrocksopp et al.l l2006f ) . It is 
very difficult to make a quantitative comparison of the 
winds seen by Chandra and XMM-Newton because the 
Doppler widths of the lines are so different: a > 3000 km 
s _1 in the XMM-Newton spectra and a « 600 km s _1 in 
the Chandra observations. This difference may influence 
the apparent column dens ity and ionization param eters. 
Nevertheless, Figure 5 of iDi'az Trigo et al.l pOOl indi- 
cates that between March 18 and March 27, the absorber 
began t o resemble the forest of lines found bv lMiller et al.l 
(2006aJ). In this context, it is also an interesting coinci- 
dence that the blucshift of the Fe xxvi line in observation 
1 is the same as the blueshift of the slow MHD wind in 
observation 2. 

Thus it is possible that a precursor to this MHD wind 
was present during the first Chandra observation. To 
explore this possibility, we expand on our XSTAR models 
from Section [3.21 In what follows, we consider a number 
of scenarios for what the properties of the precursor wind 
could have been during observation 1, or how they might 
have evolved during the outburst. To be clear, when we 
refer to a quantity as constant, we mean constant in time. 

4.2.1. Wind Density is Constant 

We begin with the assumption that the density of the 
MHD wind is constant during the outburst (n = 10 15 
cm~ 3 ), as in the xstar model shown in Figure [31 As 
discussed in Section [3.2[ this mod el is dynamically iden - 
tical to the wind in observation 2 (jKallman et al.ll2009f) . 
but features an ionization balance appropriate for obser- 
vation 1. In Section 13.21 we used this model in a sim- 
ple test to show that the wind in observation 1 cannot 
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have the same column density and ionization parame- 
ter as it does in observation 2 (Figure |3j). Here, we al- 
low iVu and £ to be free parameters and actually fit a 
warmabs model to the Chandra spectrum of observation 
1. Using a 128 x 128 logarithmically-spaced grid cov- 
ering £ = 10 3 - 10 5 and iV H = 10 22 - 10 24 , we create 
a confidence map to search for the model that best fits 
this spectrum. The absence of lower-ionization lines and 
Fe XXVI Ly/3 absorption provides upper limits on the col- 
umn density of the model and lower limits on the ioniza- 
tion parameter, respectively. We allow the MHD precur- 
sor wind to contribute as much as possible to the Fe XXVI 
absorption line and assume the rest comes from the ther- 
mal wind. If the column density is too low, the predicted 
iron line will be too weak. Furthermore, because the col- 
umn density is constrained from above, the ionization 
parameter cannot be arbitrarily high and still contribute 
to the line. Thus the Fexxvi line provides lower lim- 
its on the column density and a upper bounds on the 
ionization parameter. For n — 10 15 cm -3 , we find an 
equivalent hydrogen column density Nn = (9 ± 4) X 10 22 
cm -2 and an ionization parameter log£ = 4.2^q 2 ; the 
errors are 99% confidence limits for a single parameter. 
This result implies that the column density of the MHD 
wind increased by an order of magnitude between the 
first and second Chandra observations. 

But the column density cannot be the only parameter 
that changes between these two observations, since 



N H = n AR, 



(3) 



where AR is the (integrated) radial extent of the ab- 
sorber. In general, in order for the column density to 
change, there must be variation in the gas density, the 
extent of the wind, or both. That is, in order to explain 
the differences in the Chandra lines, we require a signif- 
icant change in at least one of these properties. In the 
model described in this section, the density is constant, 
which in turn implies from Equation [3] that AR must 
have been an order of magnitude smaller during observa- 
tion 1 than it was 20 days later. Although this scenario 
is not impossible, it is unclear why the wind would have 
the same density but arise from a much smaller region 
of the disk or have a much smaller filling factor (i.e. why 
it would be so much clumpier). In other words, we can- 
not rule it out, but we are hesitant to accept a model 
in which the column density changes solely because of 
changes in AR. 

4.2.2. AR is Constant 

At the other extreme is the case where AR is constant, 
and the change in the wind column density between the 
two Chandra observations is due solely to a change in 
gas density. In this case, it is less straightforward to find 
the density of the wind during observation 1, because 
the ionization balance and resulting fits with warmabs 
are sensitive to the density, meaning that an iterative 
process is required. Ideally, to be fully self-consistent, we 
need the density n M used to create the warmabs model to 
be the same as the density uar implied by Equation [3l 
i.e.: 



AT] 



H.l 



N- 




lo Slo( n w Cm3 ) 

Fig. 6. — Densities (top) and radii (bottom) implied under var- 
ious assumptions by our grid of warmabs models (at densities n w ) 
for the MHD precursor wind in Chandra observation 1. The open 
symbols correspo nd to the case where AR is constant during the 
outburst (Section l4. 2.21 1. while the filled symb ols correspond to the 
case where AR/ R is constant (Section 14.2.3^ ■ The solid lines are 
polynomial interpolations of these data points; the dashed line in 
the top panel represents n = n„. The intersection of an interpola- 
tion and n = n„ is a self-consistent solution for the density of the 
wind in observation 1. These self-consistent solutions are shown 
as stars. In the bottom panel, R = 0.1-Rc i s shown so that the 
implied radii may be compared to the smallest launch radius for 
the thermal component. See Section [4] for details. 

Subscripts 1 and 2 refer to Chandra observations 1 and 
2, respectively. To self-consistently determine the den- 
sity for the scenario where AR is constant, we generate 
an additional set of warmabs models, which we fit to the 
Chandra spectrum of observation 1. Again, we use the 
sa me covering fac t or, tu rbulent line width, and blueshift 
as lKallman et al.1 (|2009[ ). but we build models with den- 
sities of log(n w cm 3 ) = 11, 12, 13, 14, and 15. For each of 
these models, we perform the same grid search and con- 
fidence maps described in Section 14.2.11 the results for 
A% and £ are tabulated in Table Q] Typically, our best 
fit for any given model contributes only ^3 — 4 eV to 
the equivalent width of the Fe xxvi absorption line; the 
rest we attribute to the thermal component. 

As in the preceding subsection, all of these models re- 
quire that the column density of the MHD wind in ob- 
servation 1 be more than an order of magnitude lower 
than it was in observation 2. Since we are assuming here 
that AR is constant, the density in observation 1 must 
be lower by the same factor. The implied densities uar 
are calculated from Equation 21 listed in Table [TJ and 
plotted as open circles in the top panel of Figure [(2 the 
bottom panel of this figure shows the corresponding radii 
from Equation [T] Interpolating over our (n H , uar) grid 
with a fourth-order polynomial, we solve for the density 



at which uar = n v . We find uar 



10 



13.7 



(in- 



AR 



n 2 



H.l 



(4) 



dicated with an open star in Figure 15]). At this density, 
Equation's satisfied, so that the apparent change in the 
column density between the two Chandra observations is 
determined only by the change in density. Under these 
conditions, our constraints on the ionization parameter 
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of the precursor MHD wind in observation 1 imply a dis- 
tance from the black hole of Rar ~ 10 101 cm. This is 
outside 0.1i?c, i.e. it is consistent with the location of 
the thermal wind. 

4.2.3. AR/R is Constant 

Finally, we consider an intermediate scenario in which 
the changes in the column density are dictated by simul- 
taneous changes in n and AR. If we suppose that AR/R 
is constant, i.e. the fractional extent of the wind does not 
change during the outburst, it is easy to show that the 
implied density n A R/R is given by 



i w = "-AR/R = n 2 



Li 



A, 



H.l 



H.2 



(5) 



Again, we are primarily interested in the self-consistent 
solution, where the column densities measured from the 
data with warmabs and the density input to XSTAR satisfy 
Equation [5j so that the fractional extent of the wind is 
preserved from one Chandra observation to the next. We 
employ the same iterative approach and warmabs models 
from Section [4.2.21 for each model, we calculate tiar/r 
using our measurements of Ah, £, and L (see Table Q] 
and the filled circles in Figure |6|). The assumption that 
AR/R is constant leads to smaller densities and larger 
radii than the assumption that n or AR is constant. To 
find the self-consistent solution, we interpolate over our 
grid (filled circles in Figure ^ as before, and we find 
n A R/R ~ 10 12 - 4 cm" 3 and R AR/R « 10 10 7 cm (filled 
stars in Figure [6]). Once again, this location is also con- 
sistent with the location of the thermal wind. 

To summarize the results of this section briefly, if we 
suppose that a magnetically-driven wind makes a con- 
tribution to the FeXXVl absorption line in observation 
1, then photoionization models require that this MHD 
wind have a lower column density than in observation 2. 
Considering self-consistent geometrical interpretations of 
this difference, we find it probable that wind has lower 
density and originates farther out in the first observation; 
in many cases its location is consistent with that of the 
thermally-driven wind. Thus, in this scenario, the wind 
found in the first observation may actually be a hybrid 
MHD/thermal wind. 

5. SUMMARY AND CONCLUSIONS 

In this paper, we have explored the puzzling variability 
of the accretion disk wind of the black hole GRO J1655- 
40, as seen by Chandra during the 2005 outburst. The 
first Chandra HETGS observation, taken in March of 
that year as the source was beginning a transition out 
of the hard state, revealed a single Fexxvi absorption 
line at 7 keV; in a much softer state twe nty days later, 
Chan dra saw a forest of hundreds of lines ([Kallman et al.l 
2009), indicative of one of the densest, highest-column 
winds observed in an X-ray binary to date (so dense and 
so highly ionized that magnetic proce sses must have been 
the p rimary launching mechanism; iMiller et al.l [2006a, 
[2008h . 

But how are we to interpret the relatively sudden ap- 
pearance of so many lines? Is it possible that state- 
dependent changes in the ionizing spectrum dominate 
changes in the optical depth in the wind, or is a strongly 



variable outflow required to account for the variable 
lines? In Section 13.21 we showed that the wind could 
not possibly have been the same in both Chandra obser- 
vations, or else both spectra would have shown similar 
numbers of strong absorption lines. In fact, not only is 
the broadband spectrum in observa tion 1 incapable of 
over-ionizing the wind discovered by [M06L it should ac- 
tually be better at driving a thermal wind off the disk! 

Armed with the certainty that the wind must have 
changed during the outburst, we considered in Section 
[4] some of the ways that it might have evolved from one 
Chandra observation to the next. This was not meant to 
be an exhaustive study of every conceivable variant of the 
geometry of the wind, but a focused expl orat ion of some 
of the cleanest explanations. In Section 14. 11 we consid- 
ered the possibility that there might be two physically- 
distin ct outflows in GRO J1655-40 (see iKallman et al.1 
2009 for dynamical evidence for this conclusion) , one per- 
sistent Compton-heated wind and one transient MHD 
wind. In Section 14.21 we considered several alternate 
scenarios in which both the thermal component and the 
MHD component were present in both observations; we 
modeled the MHD component with a grid of xstar mod- 
els to constrain its properties. In these scenarios, the 
absence of so many lines during the first Chandra obser- 
vation indicates a significantly smaller column density for 
the precursor MHD wind, which could be due to changes 
in its radial extent, thickness, or filling factor, its den- 
sity, or both. In the latter two cases, the location of the 
MHD component is consistent with the location of the 
thermal component (implied by the measured ionization 
parameter and Compton temperature, respectively). 

Although there is no compelling empirical evidence 
highlighting any of these scenarios as preferable, there 
is tantalizing evidenc e from contemporaneo us XMM- 
Newton observations (jDiaz Trigo et al.l 120071 ) that the 
single-line wind detected in the first Chandra observation 
might evol ve smoothly into th e rich, dense outflow dis- 
covered bv lMiller et al.l (|2006a[) . with a se cond, more ion- 
ized a nd faster component mentioned by Kallman et al.1 
(2009). Following this line of reasoning, we suggest that 
the first Chandra observation might be probing a hy- 
brid MHD / thermal wind from the outer disk that evolves 
over the next three weeks into two distinct flows with 
very different properties. Future theoretical work on 
the expected behavior of disk winds in outburst may 
shed more light on these results, and whether thermal 
wind models can ac tually produce the observed lines (see 
iLuketic et al.l [2010). For the moment, it is remarkable 
that such insights, even if tentative, can be gleaned from 
a spectrum with only one line. But we have the benefit 
of a great deal of context from other observations during 
the outburst, as well as high-quality spectra from RXTE. 

Regardless of the precise physical explanation, it 
is abundantly clear that the differences between the 
two Chandra spectra cannot be explained by differ- 
ences in the number of ionizing photons. Instead, the 
changes in the wind are very likely linked in a fun- 
dam ental way to the state of the accretion disk (see 
alsolLee et all 120021: iSchulz fc Braiidt|l2002t IMiller et all 
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Neilsen et ail 120111 120121) . It would be particularly in- 



teresting if the changes in the magnetically-driven wind 
could be tied to changes in the magnetic field configu- 
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ration in different states, but different accretion rates, 
radiative efficienc ies, and scale heights in the inner ac- 
cretion disk (e.g. iNeilsen et aHl2012t iPonti et all 120121) 
may also be important. In a forthcoming paper, we will 
perform a comprehensive study of the radio, infrared, 
and X-ray emission in order to explore the relationship 
between the accretion state and the wind in much more 
detail. 
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